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The recently discovered, extremely proton-rich nuclide '*Mg exhibits ground-state decay via two sequential 
two-proton (2p) emissions through the intermediate nucleus, 1°Ne. This study investigates the structure and the 
initial 2p decay mechanism of !*Mg by examining the density and correlations of the valence protons using a 
three-body Gamow-coupled-channel method. The results show that the ground state of '*Mg is significantly 
influenced by the continuum, resulting in a significant s-wave component. However, based on the current frame- 
work, this does not lead to a significant deviation in mirror symmetry in either the structure or spectroscopy of 
the '*Mg-'8C pair. Additionally, the time evolution analysis of the ‘Mg ground state suggests a simultaneous 
2p emission during the first step of decay. The observed nucleon-nucleon correlations align with those of the 
light-mass 2p emitters, indicating a consistent decay behavior within this nuclear region. 
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I. Introduction 


In recent decades, remarkable advancements have been 
made in extending the boundaries of nuclear stability, lead- 
ing to the discovery of numerous dripline nuclides. These nu- 
clides, which are characterized by a significant imbalance in 
their proton—neutron ratios, exhibit exotic modes of radioac- 
tivity and are significantly influenced by the continuum ef- 
fect. Notably, two-proton (2p) emission, initially proposed 
by Goldansky in the early 1960s [1, 2], has only been ob- 
served in recent years [3, 4], making it a particularly rare 
form of radioactivity. This decay process is essentially a 
three-body quantum-tunneling phenomenon occurring in ex- 
tremely proton-rich nuclei. It indicates the complex interplay 
between the daughter nucleus and valence protons [5-9]. The 
decay mechanisms vary across different 2p emitters [10-14], 
encompassing diproton decay and large-angle emission [15- 
18], depending on the internal structure of the nuclei and the 
properties of adjacent nuclei [19-23]. This diversity in de- 
cay processes and the potential to investigate them via proton 
emission correlation measurements have garnered significant 
interest and serve as crucial testing grounds for both nuclear 
physics and quantum theory. 

Investigations beyond the proton dripline have revealed nu- 
clei that are unbound with respect to new decay channels, 
including multi-proton emitters. This highlights the signifi- 
cance of such exotic decay modes as the precursor decay en- 
ergy increases. Notable examples include nuclei that emit 
three protons (3p), such as 7B [24], !”Na [25], 3!K [26], and 
13F [27], as well as ground-state emitters of four protons (4p), 
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like 8C [28] and !8Mg [29], and a five-proton (5p) emitter, 
°N [30]. Although the decay characteristics of many of these 
systems have not been investigated, analyses show that their 
decay may proceed in multiple steps, shedding excess protons 
through sequential 1p or 2p emission. "Mg, which has been 
recently discovered [29], decays via sequential 2p—2p emis- 
sion of two pp pairs. Therefore, studies investigating whether 
the first-step 2p decay of !8Mg resembles that of other 2p 
emitters and how the inner structure impacts this process are 
required. 


Addressing the complexities of these phenomena requires 
sophisticated theoretical frameworks; however, comprehen- 
sive models that can simultaneously address the nuances of 2p 
and multi-proton emissions are still under development. Var- 
ious approaches [31—34]have been employed, including the 
configuration interaction and density functional methods for 
structure analysis [35—40], and few-body methods for probing 
asymptotic proton correlations [41—45]. Innovations such as 
the shell model embedded in the continuum and the Gamow 
shell model integrate continuum effects into the configura- 
tion interaction framework, whereas hybrid methods leverage 
spectroscopic factors and few-body partial decay widths. Al- 
though the spectra and isospin-symmetry breaking between 
18Mg and mirror symmetry ‘SC have been extensively in- 
vestigated [46-49], the properties of the SMg-!8C pair, es- 
pecially the structures and time-dependent nucleon-nucleon 
correlations, have not yet been elucidated. This study em- 
ployed the Gamow coupled-channel (GCC) method, a three- 
body approach that incorporates continuum effects, to inves- 
tigate the decay dynamics and asymptotic correlations of 2p 
emitters, with a particular focus on the first-step 2p emission 
of '8Mg and the configurations of its valence protons. 


The remainder of this paper is organized as follows. Sec- 
tion II describes the models and Hamiltonian. In particu- 
lar, it presents the framework of the GCC method and time- 
dependent approach. The structural and decay properties of 
18Mg are discussed in Section III. Finally, the summary is 
presented in Section IV. 


Il. Theoretical Model 
A. Gamow coupled-channel approach(GCC) 


In this study, we used a GCC approach [50, 51], which is 
a three-body model describing !8Mg in terms of two valence 
protons and a core (‘°Ne). To describe the three-body asymp- 
totic behavior, the coordinates of !8Mg can be expressed in 
Jacobi (relative) coordinates: 


cogs hs Aba gas r. (1) 
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where r; denotes the position vector of the ith cluster; 
i1 = pı, i2 = po, and ig —!°Ne for the T-coordinate and 
i, =l6Ne, i2 = pı, and iz = po for the Y-coordinate, as 
shown in Fig. 1(a,b). A; is the mass number of the ith clus- 
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masses associated with x and y, respectively. These parame- 
ters can be used to describe the hyperradius p = y£? + y?, 
which is transformation-invariant among different sets of Ja- 
cobi coordinates. 

As the experimental measurements are performed in the 
momentum space, the relative momentum is defined as fol- 
lows: 
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where k; denotes the linear momentum of the ith cluster; 0; 
and 0; are the opening angles of (kx, ky) in the Jacobi-T and 


Jacobi-Y coordinates, respectively, as shown in Fig. 1(c), 
24.2 


respectively. Epp = The represents the kinetic energy of 
the relative motion of the two emitted protons, and Ecore—p 
signifies that of the !°Ne-p pair. 

The total wave function in the Jacobi coordinates of the 
mother nucleus, which couples with the p——p pair and ‘°Ne, 


can be expressed as 
wt = S [o’" Q prene |r, (3) 


where 7°" denotes the wave function of the core 1ŝNe, and 
QIe™ represents the wave function of the valence nucleons, 
which can be expressed as 
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where y (0) reprsents the hyperspherical harmonic of 
the hyper-angle. The hyperradial part is expanded in 
the Berggren ensemble, a complete basis in the complex- 
momentum plane including bound, decaying, and scattering 
states [50, 52, 53]. K is the hyperspherical quantum number, 
and y = {51, $2, S12, S, lx, ly, L, Jp, Je} contains the other 
numbers. Using the Berggren basis, the inner and asymptotic 
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Fig. 1. Schematic of three-body '°Ne + p +p system in Jacobi-T (a) 
and Jacobi-Y (b) coordinates. Momentum scheme (c) of three-body 
system. 


regions of the Schrédinger equation can be addressed on the 
same footing, providing a natural connection between the nu- 
clear structure and decay aspects of the problem. 

The !©Ne + p + p Hamiltonian of GCC can be expressed as 
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where p? /(2m,;) represents the kinetic operators of each par- 
ticle, Vij (rij) denotes the interaction between clusters 7 and 
J, H, is the core Hamiltonian given by the excitation energies 
of the core, and OE ee denotes the center-of-mass term. In 
Jacobi coordinates, the center of mass is automatically elimi- 
nated. In this study, the proton-core interaction Vpc, including 
the central, spin—orbit, and Coulomb terms, is approximated 
by a Woods-Saxon (WS) average potential[54—57]. 

To address the antisymmetrization between core and va- 
lence protons, a supersymmetric transformation method [58- 


60]was employed. This method introduces an auxiliary repul- 
sive “Pauli core” in the original core-p interaction to eliminate 
Pauli-forbidden states. 

Consequently, the acquired eigenstate possesses a com- 
plex energy denoted by EF = E — il /2, where I repre- 
sents the decay width. To elucidate the dynamics and asymp- 
totic correlations inherent in the 2p decay, the wave func- 
tion WE derived within the complex framework, should 
be decomposed into real-momentum scattering states by ap- 
plying the Fourier—Bessel series expansion. This approach 
is instrumental in reinstating the Hermitian characteristics of 
the Hamiltonian matrix and ensuring the conservation of total 
density. The initial wave packet, W47 (t = 0), thus generated, 
can evolve over extended durations and distances within a 
time-dependent (TD) framework, as shown in Refs. [61, 62]. 


B. Hamiltonian and parameters 


Magnesium isotopes often exhibit significant deformation. 
However, due to the lack of definitive experimental data on 
their deformation, the deformation factor is adopted from 
their mirror nucleus, 18C, which exhibits a quadrupole de- 
formation, 62, of approximately 0.3, as reported in [63]. The 
ground and 2* states of the deformed core were considered 
through non-adiabatic coupling with the valence nucleons, 
and the corresponding energies were obtained from experi- 
mental data [63]. A previous study [38] reported significant 
mirror symmetry breaking in the ‘°Ne-!°C mirror pair due 
to different configurations [64, 65]. However, the excitation 
energies of their first 2+ states are remarkably similar, with 
16Ne at 1.77(0) MeV and 16C at 1.76(6) MeV. Consequently, 
the breaking of mirror symmetry within the current frame- 
work cannot easily be described. 

In a three-body system (1°Ne + p + p), the interaction be- 
tween the two protons is described by the finite-range Min- 
nesota force, adhering to the original parameters specified 
in Ref. [66]. This interaction is further supplemented by a 
two-body Coulomb force applied to protons. The core-proton 
(core — p) effective interaction was modeled using a WS po- 
tential (with a spin-orbit term) and one-body Coulomb inter- 
action. The parameters of the WS potential, except for the 
potential depth Vp and- the spin-orbit strength V,,., were re- 
tained as default values, as detailed in Ref. [67]. The spin- 
orbit strength was set to 24 MeV. Additionally, the depth Vo 
was finely adjusted to accurately replicate the experimental 
ground state decay energy, Q2p(0" g.s.) = 3.44 MeV, and the 
energy of the first 2+ state, Ez (27) = 1.84 MeV, of !8Mg, as 
discussed in Ref. [29]. 

The three-body configurations in Jacobi coordinates are la- 
beled by quantum numbers (K, Lx, ly, S), where ly repre- 
sents the orbital angular momentum of the proton (neutron) 
pair with respect to their center of mass, £, is the orbital an- 
gular momentum of the pair with respect to the core, and S 
denotes the total intrinsic spin of the emitted nucleons. Three- 
body calculations were performed in the model space defined 
by max (lz, ly) < 8 with a maximal hyperspherical quantum 
number Kmax = 20. For the hyperradial part, we used the 
Berggren basis for the K < 3 channels and the harmonic os- 


cillator basis with an oscillator length of 1.75 fm and Nmaz 
= 20 for the remaining channels. For the GCC calculation of 
the initial state, the complex momentum contour defining the 
Berggren basis is given by the path k = 0 > 0.3 — 0.li — 
0.4 — 0.052 — 0.5 — 0.037 + 0.8 4 1.2 => 2 —> 4—6 
(all in fm~). For time-dependent evolution, the inner part (< 
15 fm) of the initial state was expanded and propagated with 
a real-momentum contour, which followed k = 0 > 0.1 > 
0.13 > 0.15 — 0.16 — 0.17 > 0.185 > 0.2 — 0.225 > 
0.24 > 0.25 > 0.65 > 0.8 31.2525 4 > 6 —> 8 all 
in fm~'). Each segment was discretized into 30 scattering 
states. In practice, interactions within a sphere of radius 250 
fm were considered. However, because the wave function is 
defined in the momentum space and evolves from a highly lo- 
calized initial wave packet, this cutoff has no practical effect 
on the investigated physical observables. 


Ill. Results and Discussions 
A. Structure information of '*Mg and ‘°C 


A highly proton-rich nucleus, ‘Mg, located far from the 
B-stability line, has been recently identified as a 4p emitter 
[29]. The 4p decay proceeds through a two-step 2p decay 
via the intermediate ground state of '©Ne. To elucidate the 
internal structure of !°Mg and its influence on decay proper- 
ties, we begin by examining the spectrum and valence-proton 
configurations of "Mg and comparing them with its mirror 
partner, 18C. The mirror pair, due to isospin symmetry, typ- 
ically exhibits similar properties across many aspects. How- 
ever, Coulomb interactions introduce discrepancies in their 
thresholds. Consequently, !8Mg is particle-unbound, with a 
2p separation energy of S2p(18Mg) = —3.44 MeV, whereas 
18C is bound to a 2n separation energy of S2,,(18C) = 4.92 
MeV, as determined experimentally. 

The calculated spectra of Mg and 18C align qualitatively, 
if not quantitatively, with the experimental data. Specifically, 
the excitation energy of the first 2+ state is slightly under- 
estimated for "Mg. Meanwhile, although only one 2* state 
is observed experimentally, our calculation predicts two low- 
lying 2+ states, similar to the situation in its mirror system 
18C, The ground state of !85Mg is less unbound than those 
of !Na and '®Ne. Considering the large decay width of 
17Na, the ground state of '*Mg might have a “democratic” 
decay mode, which exhibits open 1p and 2p decay channels. 
This suggests potential competition between 1p and 2p de- 
cay modes during the decay process of '?Mg. To determine 
the primary decay mechanism, we used a time-dependent ap- 
proach to analyze the corresponding decay dynamics (see the 
discussion below). 

However, due to the Coulomb interaction, 18Mg is less 
bound than its mirror partner, indicating a greater continuum 
effect. As presented in Table 1, the sd-shell effect signifi- 
cantly influences the ground-state configuration of valence 
neutrons in 18C, favoring the configurations (K,lz, ly, S) = 
(4,0,0,0) and (4,1,1,1), which primarily correspond to the d 
orbital. However, for !85Mg in its ground state, the (0,0,0,0) 
configuration (s orbital) is essential. This is because 18Mg 
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Fig. 2. The calculated spectra and decay widths (shaded areas) of '*Mg-'8C pair and related neighboring nuclei "Na and '’C. Experimental 


data were obtained from Refs. [29, 63]. 


TABLE 1. Predicted configurations (K,éz,€y,S) and the corre- 
sponding weights of Mg and ‘SC. 


J” (K, lx, ly, S) 

OR. 34.86% (4,0,0,0) 39.68% (4,0,0,0) 
19.36% (4,1,1,1) 25.82% (4,1,1,1) 
18.00% (0,0,0,0) 6.83% (4,0,2,0) 

git 13.36% (4,0,2,0) 17.28% (4,1,1,1) 
11.61% (2,0,2,0) 14.10% (4,0,0,0) 
8.29% (4,0,0,0) 13.59% (4,0,2,0) 

at 25.51% (4,2,2,0) 25.57% (4,1,1,1) 


16.41% (4,1,1,1) 
4.39% (8,4,4,0) 


13.67% (4,2,2,0) 
11.06% (4,0,0,0) 


2t 23.39% (4,0,0,0) 
15.51% (4,1,1,1) 
9.15% (2,0,2,0) 


better interacts with the continuum than 18C, rendering the 
s-wave component favorable owing to its reduced centrifu- 
gal barrier. This effect is also observed in the spectral anal- 
ysis of the odd-A neighboring nuclei ‘Na and 17C, where 
the 1/2* state in '’Na, despite being unbound, benefits from 
additional binding from the continuum, establishing it as the 
ground state. Conversely, the state sequence in the mirror nu- 
cleus 17C exhibits a different arrangement. This discrepancy 
in the energy levels between particle-bound and unbound mir- 


ror systems is a characteristic of the Thomas—Ehrman shift 
[68-70]. 

Regarding the 2* states, although the pure s-wave compo- 
nent is forbidden by the angular momentum selection rule, the 
influence of the continuum is still noticeable, causing the con- 
figuration (K’,¢,,, ly, S) = (2,0,2,0) to be pronounced in ‘Mg, 
where the K quantum number represents the three-body cen- 
trifugal barrier. This leads to similar excitation energies for 
the 2+ states of 18C and !8Mg, although different continuum 
effects are involved in these systems. 

The internal structures of 18C and !°Mg also exhibit mod- 
erate mirror symmetry breaking, which is similar to that of 
the 16C-!6Ne mirror pair. As illustrated by the density distri- 
butions of the valence nucleons shown in Fig. 3, the ground 
state of !8Mg exhibits a more expansive distribution than that 
of 18C, although both nuclei possess di-nucleon, cigar-like 
structures, and a small triangular component internally. This 
discrepancy is due to the significant s-wave component and 
unbound nature of !SMg. Meanwhile, the internal structures 
of the a states appear analogous to the ground for both “Mg 
and 18C because most components are built upon the excited 
state of the daughter nuclei. 

The properties of the related nucleus, !’Na, have also at- 
tracted considerable attention due to limited experimental 
data availability. Energy spectra have been calculated in sev- 
eral studies to predict the decay scheme of '’Na [71, 72]. 
In addition, to assess the effect of deformation on structural 
characteristics and decay properties, we computed the energy 
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Fig. 3. 2p (left) and 2n (right) density distributions in Jacobi-T 
coordinate predicted for the ground and first excited states of !8Mg 
(left) and 18C (right), respectively. 
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Fig. 4. Calculated spectrum of "Na as a function of quadrupole 
deformation 62. 


spectrum of !”Na and the 2p density distributions of the !8Mg 
ground state as functions of the quadrupole deformation pa- 
rameter (2 (refer to Figs. 4 and 5). The analysis of the spec- 
tral evolution of "Na reveals a gap between the s and d or- 
bitals with increasing deformation. This gap is also observed 
as the (K,lz, ly, S) = (4,0,0,0) and (4,1,1,1) configurations 
decreases from a spherical shape to a prolate deformation 
in the ground state of !8Mg. Consequently, this results in 
a change in the primary structure of !85Mg, as shown in Fig. 
5. Notably, considerable deformation (82 = 0.5) represents 
uniform density distribution. 


B. Decay properties of '“Mg 


The internal configurations and spectroscopic details of 
18Mg and its adjacent nuclei have been reported previously. 
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Fig. 5. Ground-state density distributions of *5Mg with (a) spherical 
and (b) 82 = 0.5 deformations. 
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Fig. 6. Time evolution of the ground state of '*Mg with a quadrupole 
deformation G2 = 0.3. Density distributions of two emitted protons 
are shown in Jacobi-Y (left) and -T (right) coordinates for four dif- 
ferent time slices. 
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Fig. 7. Asymptotic energy (a), angular (b), and two-dimensional cor- 
relations of two emitted protons from the ground state of *Mg with 
{2 deformation of 0.3 in Jacobi-T (left) and -Y (right) coordinates, 
respectively. Q2, is the total decay energy. 


An intriguing aspect to further examine is the influence of 
these properties on the decay dynamics and mechanisms. For 
this purpose, we employed a time-dependent framework to 
evolve the initial wave function across extensive spatial and 
temporal scales, as shown in Fig. 6. 


Initially (t = 0), most of the !85Mg ground-state wave func- 
tion was confined within the nuclear boundary, forming the 
discussed diproton and cigar-like structures. However, the 
density evolution within the Jacobi-T’ coordinate reveals the 
merging of these structures during the tunneling phase; a dy- 
namic distinct from the decay behavior observed in the p- 
shell nucleus Be [61] but similar to that of the sd-shell nu- 
cleus !7O [73]. This similarity is attributed to the substantial 
s-wave component and relatively large decay energy/width 
characteristic of the !8Mg ground state. Concurrently, in the 
decay process, the distance between '°Ne and a proton re- 
mains similar to the increasing separation in the (1°Ne,p)-p 
configuration, as observed in the density evolution depicted 
in the Jacobi-Y coordinate shown in Fig. 6. This observation 
suggests that, despite the energetic feasibility of single-proton 
decay for the !85Mg ground state, simultaneous two-proton 
emission is the most probable decay mode. 

The correlations between protons emitted from the ground 
state of Mg are shown in Fig. 7. Given the notably wide 
decay widths associated with light—mass systems, the energy 
or angular distributions cannot easily reveal distinct patterns 
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Fig. 8. Similar to Fig. 7 but for the $2 deformation of 0.5. 


that could differentiate between decay mechanisms. Specif- 
ically, for Mg, the energy correlation of the T-type is al- 
most uniform, which is attributed to the substantial s-wave 
component within the ground-state wave function, as previ- 
ously discussed. Meanwhile, the Y-type energy correlation 
attains a noticeable peak around Eyore—p / Qo», further sup- 
porting the likelihood of a simultaneous 2p emission, rather 
than a sequential decay process. 

Furthermore, we investigated the effect of deformation on 
the asymptotic correlations of “Mg, as shown in Figs. 7 
and 8. When the quadrupole deformation 82 changes from 0 
to 0.3, the asymptotic correlation is slightly affected. How- 
ever, the T-type energy correlation and Y-type angular cor- 
relation are significantly changed for 82 = 0.5. This is in 
accordance with the change in the internal structure/configu- 
ration as deformation increases. In particular, enhanced peaks 
are observed in the small Epp and large0; regions (see Fig. 
8), which might indicate a diproton emission, although the 
density of the valence proton inside the nucleus is uniformly 
distributed for 62 = 0.5 (see Fig.5). Given that nucleon— 
nucleon correlation can be directly measured, further experi- 
mental investigations are required. 


IV. Summary 


In this study, we investigated the structure and 2p decay 
mechanism of !85Mg using a three-body Gamow-coupled- 
channel model and compared it with its mirror partner, 18C. 


Our analysis revealed that although both systems exhib- 
ited di-nucleon and cigarlike structures in their valence nu- 
cleon configurations, moderate mirror symmetry breaking oc- 
curred in the density distributions and spectroscopy due to the 
Thomas—Ehrman effect. The influence of the continuum on 
the ground state of '*Mg induced a significant s-wave com- 
ponent. 


The structural characteristic of "Mg led to widespread 
distribution of the nucleon-nucleon correlation within the 
asymptotic region, similar to the behavior of other 2p emitters 
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